Suppression of modes in the random phase approximation 
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A general but simple method is proposed to eliminate the quantum fluctuations generated by 
selected one-body operators in the excitation spectrum of a discrete RPA Hamiltonian. This method 
provides an outstanding tool for the removal of the contaminating spurious effects originated from 
symmetry violations. It can be also applied as a mode filter for analysising RPA response functions. 

PACS numbers: 21.60.Ev, 21.60.Jz, 24.30.Cz 
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The Random Phase Approximation (RPA) P, |3| is 
a powerful standard approach to calculate microscopi- 
cally the variety of vibrational excitations and giant res- 
onances in nuclei and in other finite fermion systems like 
quantum dots or metal clusters 0, S 0|- Based upon 
a stable equilibrium configuration of independent quasi- 
particles in a mean field (MF) potential, the RPA is ac- 
counting for the quantized small amplitude oscillations 
(named below as fiuctuations) about the MF equilibrium 
point by allowing Bose-like two-quasiparticle (2qp) exci- 
tations driven by the residual interaction. There is rich 
diversity of phonon excitation modes produced in this 
way which shows up as complex resonance pattern in 
calculated RPA response functions. 

In this letter we present a new method which enables a 
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FIG. 1: Illustration of the suppression of the spurious cm. 
fluctuations in the RPA response function of the isoscalar 
dipole operator D = r^Yim for a nuclear system with equal 
proton and neutron numbers Nt^ = Nv = 20. The RPA 
Hamiltonian consists of a spherical Nilsson potential and a 
residual interaction of dipole plus octupole type. The selfcon- 
sistent value Ksc of the isoscalar dipole strength is mul- 
tiplied by an increasing factor in order to demonstrate the 
disappearance of fluctuations SR^ of the cm. coordinate R. 
In the solid curve the fluctuations of both the cm. coordinate 
R and the cm. momentum P are damped out. 
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FIG. 2: Mode suppression effects illustrated in a schematic 
RPA calculation of the Ml response function for nucleus 
^^^Gd. The RPA Hamiltonian consists of a rotating axially 
deformed Nilsson plus pairing potential and a residual interac- 
tion of spin plus quadrupole type. Solid curve: Fluctuations 
of both the isoscalar (T = 0) a.m. J„ + J„ and of the par- 
ticle numbers NTr,N„ suppressed. Dotted curve: without any 
suppression. Dashed dotted curve: Like the solid curve but 
additionally the fluctuations of the isovector (T = 1) a.m. 
part Jtt — Jv suppressed. Note, that in the latter curve the 
dramatic effect might be not realistic because of the simple 
interaction. 



forced damping of selected degrees of freedom. On the 
one side this method provides us a perfect tool for the 
elimination of spurious modes as e.g. the unwanted cen- 
ter of mass (cm.) motion as examplified in Fig. 1 or the 
fiuctuations of the angular momentum (a.m.) and the 
particle number as demonstrated in Fig. 2 which anyway 
contaminate the calculated response functions. On the 
other side the method can be used as a more general in- 
strument for probing system properties by switching off 
the fluctuations of selected excitation modes for a spec- 
tral analysis as demonstrated in Fig. 2 for a single mode. 
The underlying idea of our method is quite simple. Let 
us consider a one-body operator F which by definition 
incorporates a specific mode of 2qp excitation. Then, by 
adding the oscillator-like restoring force term kF^ to the 
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considered RPA Hamiltonian, one can shift arbitrarily 
the oscillator frequency by tuning the strength param- 
eter K. In the limit k ^ oo the restoring force gets so 
large that the RPA fluctuations of the operator F become 
completely suppressed, i.e. the operator F is frozen at 
its constant MF value. In this spirit the excitation of the 
mode F gets eliminated from the RPA excitation spec- 
trum. Hence, by choosing appropriate force terms the 
corresponding excitation modes can be systematically fil- 
tered out (cf. Fig. 2). 

Now the essential steps of our derivation will be sketched 
using well known relations P, ■ 
Let us consider a whole group of one-body operators P{r) 
and Q{r) selected for the removal of the associated modes 
(denoted below as PQ modes). Examples for such op- 
erators are the cm. coordinates R, the particle number 
operator N or the multipole operators Qim — r^Yi„i. The 
PQ operators are given by 



P[t) 



Q{r) 



Y^c^^(^ + h^ (1) 



where we restrict ourselves without loss of generality to 
hermitean coordinate-like and momentum-like operators 
with any real coefficients p^''^ and g'-''-*. In eqs.(l), bi and 
denote the basic 2qp excitation operators treated in 
RPA as bosons obeying approximately the quasi boson 
approximation (QBA) [6^,6^"] « Sij In practice, 

the explicit bosonic form, eq.(l) follows from the original 
fermion representation of the Ps or Qs by a qp transfor- 
mation and using subsequently the QBA. Possible non- 
zero expectation values of the quasiparticle operator like 
the average particle number can be omitted since they 
are irrelevant constants. 

Let be H{b, 6+) the bosonic Hamiltonian of the system 
under study derived analogously by applying the QBA. 
As mentioned above the key of our removal procedure is 
to add an oscillator-like force term to H{b, 6+) in which 
the sum of the squared PQ operators is collected: 



H'{k) = H{b, b+) + KY,[P{r? + Q{r) 



(2) 



The RPA excitation modes s and their energies E'^ are 
constructed by solving the equations of motion as usual 
but here for the supplemented Hamiltonian H': 



[H'in),0+is)]^EiO+{s), 
0+(,s)^5][X,(s)6+-r,(s)6,]. 



(3) 



The operators 0^(s) denote the familiar excitation oper- 
ators of the RPA eigenmodes and the amplitudes {X, Y) 
are built from the eigenvectors which diagonalize the 
well-known RPA matrix to H'{k). It is irrelevant for 
our derivations that in general not all modes s do have 
normalizeable amplitudes (cf.jlj). 



The tunable PQ force term in the Hamiltonian H', eq.(2) 
serves for shifting the selected modes away from the en- 
ergy region of the usual RPA excitations in order to 
achieve the intended mode suppression. In this spirit 
our method extends the known idea of pushing out un- 
wanted excitation modes from the physically interesting 
spectrum as it was proposed, e.g., for the elimination 
of the cm. motion in order to restore the Galilean in- 
variance 0|. Similar requirements were formulated with 
respect to the rotational spurious motion fl^. However, 
these options have not been pursued consequently in the 
past. 

The common strength parameter k in eq.(2) is used as a 
scale variable to achieve the crucial spectral shift of the 
PQ modes. The scaled Hamiltonian H" given by 

i7»^^^^i^+5:[PMVQ(rf] (4) 



has obviously the same solutions 0^(s) as H' to the 
scaled eigenvalues i?" = E'^/ k. In the limit k ^ oo the 
term H{b,b^) / k vanishes asymptotically and the spec- 
trum of H" is determined alone by the remaining supple- 
mentary term. (Note, that this statement and its conse- 
quences below would be the same for the limit n — > —oo.) 
The RPA eigenvalues E'^ can be either exactly zero or 
have any real or complex values. For a sufficiently large 
value K = Umax < CO the part of the eigenvalue spec- 
trum related to the scaled term H(b,b'^)/ k will come 
close to zero whereas the dominating non-scaled sup- 
plementary term ^j.[P{r)'^ + Q{r)'^] has possibly exact 
zero eigenvalues (see below) or other eigenvalues but 
those certainly converging to a finite distance to the 
zero point. Accordingly, the corresponding spectrum of 
the non-scaled Hamiltonian H' , eq.(2) decomposes after 
rescaling Ei. — kE'J into three parts: (i) the exactly zero 
eigenvalues, (ii) eigenvalues shifted to infinity both re- 
lated to the supplementary term in eq.(2) and (iii) the 
"normal" vibrational eigenvalue spectrum related mainly 
with the original part H{b, b'^) to be cleaned up from se- 
lected predefined PQ modes. Therefore, by including a 
supplementary term and scaling we succesfuUy separate 
in the asymptotic region k = Kmax < oo the normal 
phonon spectrum from the part of the PQ modes. 
The resulting separation of the energy spectrum for large 
enough values k translates in a partitioning of the RPA 
solutions. We denote the asymptoptic zero and infinity 
energy solutions as 0(s = a) = Oa and 0'^{s = a) = il+ 
and the normal vibrational solutions belonging to the fi- 
nite energies as 0(s = n) = a;„ and 0~^{s = n) — a;+. 
Since by construction the sets fla, ^-nd Wn, i repec- 
tively, belong to different energy regions these modes are 
orthogonal, i.e., they satisfy the orthogonality relations 
OJ: 



0. 



(5) 
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Defining |0) to be the normalized vacuum state to the 
normal solutions u}~^{n), the vibrational states \n) are 
created by the phonon excitation 

\n)^u;+\0). (6) 

Hence by using eq.(5) one obtains 

mt\n) = main) = (7) 

for all a and n, respectively. 

We recall that in the asymptotic region k — > oo the 
set ila,^a can be considered as a RPA solution to 
the scaled Hamiltonian H" = J2r[P('^f + eq.(4) 
such that one can construct these asymptotic solutions of 
[H", = Ea^a directly in terms of the PQ operators 

^Y.^Xr{a)P{r)+Yria)Q{r)] (8) 

r 

instead of solving the full equations of motion (2). Re- 
lying on the existence of a closure relation 1] the above 
equation can be considered as formally inverted to get 
the expansions 

a 
a 

Hence, by using eqs.(7,9) we arrive at the crucial relations 

(0|F(r)|n) (0|Q(r)|n) =0 (10) 

in the asymptotic region n —^ oo . These relations express 
that all vibrational transition amplitudes for any selected 
PQ operator vanish, i.e. all those operators cannot create 
vibrations from the phonon vacuum. Thus, it is proven 
that the inclusion of a supplementary interaction term 
into the boson Hamiltonian enables one to eliminate the 
fluctuations caused by all the operators P(r) and Q{r) 
from all the normal solutions tj^. In other words, the 
response function of any selected PQ operator is identi- 
cally zero in the physical spectral region. 
Defining the sum rule Sp of a linear boson operator 
F (i.e. F is a 2qp excitation operator) as Sp = 
E„l(0|^l?i)P = (0|-F^|0) = SF^ then eq.(lO) can be 
written in compact form as 

i.e. the complete disappearance of vacuum fluctuations 
in the asymptotic region n oo for all P and Q. The 
above result (11) can be interpreted as a total blocking 
of RPA fluctuations such that the momenta P and coor- 
dinates Q get fixed. 

For illustration we present two simple but typical exam- 
ples for the supplementary Hamiltonian (2). The first 



example is the forced damping of the fluctuation for a 
single operator P (or Q). Then, in the asymptotic re- 
gion holds simply H" = P^ . One immediately real- 
izes that the corresponding eigenmode is the "zero vec- 
tor" n ^ n+ ^ P with the eigenvalues E" = E = 
since [H",P] = 0. The same applies also to a se- 
ries of P— operators (or analogously Q— operators) since 
[P{r),P{r')] = [Q{r),Q{r')] = are mutually commut- 
ing. One gets a series of zero eigenvectors identical with 
the Ps or alternatively Qs, respectively. 
Choosing symmetry operators like the cm. momentum, 
the a.m. operator or the particle number, then this mode 
suppression accomplishes the successful restoration of the 
corresponding symmetries on RPA level. For instance, 
when a statical pair fleld for protons and neutrons are 
included in the MF, the fluctuations of the proton and 
neutron particle numbers can be removed by using the 
supplementary term {N^ + N^). 

The second example concerns a conjugated operator pair 
(P, Q). Defining Q — ~Q/{2 J2rPr'lr) yields the commu- 
tator [P, Q] = —i which leads to an oscillator-like supple- 
mentary Hamiltonian H" = 1/2 (P^ + Q^)- Hence, the 
asymptotic eigensolution to [H" , $7+] — E"rt~^ is explic- 
itly given by E" = 1 and n+ = l/\/2(Q + iP) which 
yields P = (f]+ - n)/V2 and Q = (S1+ + f7)/\/2 and 
E = kE" indeed shifted to infinity for k ^ oo. Taking 
e.g. for the operator pair the cm. momentum and coor- 
dinates then eqs.(ll) realize the total removal of the spu- 
rious cm. fiuctuations i.e. the cm. coordinate is fixed 
in space and the system gets held at rest (cf. Fig. 1). 
Now we consider the possible numerical problems en- 
countered in the RPA calculations when performing the 
above described mode removal in practice. Firstly, it 
should be mentioned that the mode suppression (11) by 
an asymptotic supplementary term is working for any 
finite dimensional RPA Hamiltonian independently of 
whether the inherent mean field part is selfconsistently 
calculated or not. Further, it does not matter whether 
the RPA interaction consists of a sum of factorized terms 
like the well known multipole expansion or it is non- 
factorized like the nuclear Skyrme interaction. The ac- 
tual form of the interaction only determines whether the 
eigenvalue RPA matrix equation can be performed by the 
more convenient response function method |l2l| or by a 
more involved matrix diagonalization. 
With respect to the actual performance of the transition 
K —^ Umax < oo the Critical question arises whether nu- 
merically stable results can be expected or not. The an- 
swer is clearly positive for a boson Hamiltonian H(b, b^) 
with a factorized interaction part. In this case all 
strength parameters enter inversely in the response func- 
tion matrix (cf. T^l). Hence, the asymptotic limit of 1/k 
gives a well defined zero which is causing no problem at 
all. The same should be true also for a non-factorized 
interaction no ting the latter can be expanded into a fac- 
torized series hj. In fact, we studied the convergence 
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properties of a straightforward RPA matrix diagonal- 
ization for random RPA matrices up to the dimension 
n ~ 1000 supplemented with a series of non-commuting 
random PQ terms. It turned out that the stepwise vari- 
ation of up to the considered value K„iax = 10^ leads to a 
rapidc convergence of the vanishing fluctuations (11) to- 
gether with converging RPA eigenvalues E'^ which points 
to a good numerical performance. How large the value 
Kinax should bc actually taken can be easily controlled 
by checking the descending size of the zero point fluctu- 
ations 6P^ and SQ"^. 

We are now turning to possible physical applications and 
improvements seen so far for the nuclear structure inves- 
tigations done with RPA. 

The mode suppression should be included regularly for 
the removal of spurious modes connected with violations 
of conservation rules as for instance the Galilean invari- 
ance, a.m. conservation due to a deformed and possibly 
rotating mean field, particle number conservation broken 
by a static pair potential etc. to avoid the corresponding 
fluctuations. We mention that our method is exact as 
compared to the approximate symmetry restoration \vi\ 
which in general can not satisfy the condition (11) for a 
complete suppression of the spurious motion. 
There are immediately some useful consequences of our 
method. For instance, the effective charges and Cj/ of 
the the electric dipole operator usually introduced by a 
transformation to the cm. system become obsolete 
as well as the cm. correction Qsm = f^Ysm — W^im of 
the octupole operator discussed in is not needed (cf. 
Fig. 2) when the cm. mode R is suppressed. 
For a rotating deformed system the mode suppression en- 
ables one to define cleanly a body fixed or alternatively 
an a.m. fixed reference system. The body fixed system 
is singled out by the requirement SQ2m=±i — ^Jq = 
which accomplishes that both the principal axes system 
as well as the length \ J\ of the a.m. are conserved in 
the RPA. Alternatively, the damping out of the fluctu- 
ations of the a.m. vector through the condition 6J'^=Q 
determines the lab reference by flxing length and direc- 
tion of the a.m. vector J but still leaving freedom for 
possible shape and orientation fluctuations of the rotat- 
ing deformed nucleus. The latter option could be an im- 
portant ingredient for improving the RPA description of 
wobbling excitations found in triaxially superdeformed 
rotating nuclei 13]. 

Another important new aspect is that a simplifled version 
of the RPA can be utilized where the MF part is adopted 
from the extremely successful semimicroscopic models 
that apply the Strutinsky shell correction method 0. 
The bosonic interaction part might then be added as a 



multipole-multipole interaction with adjustable strength 
parameters. The suppression of spurious modes works 
also for this case and when properly treated this simpli- 
fied RPA satisfies the stability criterion j^. 
The mode suppression enables one to perform also a RPA 
calculation based on a 2qp excited state instead on the 
MF ground state. Usually, this interesting new possibil- 
ity is excluded since excited MF configurations violate 
the RPA stability criterion However, the stability 
can be recovered by a force term that is blocking explic- 
itly the deexcitation to the ground state configuration. 
Finally, we like to emphasize that the proposed mode 
suppression causes almost no additional effort compared 
to a RPA calculation without supplementary terms. 
In summary, the restoration of symmetry violations is 
one field of applications of the proposed mode suppres- 
sion. The possible application as a spectral filter offers as 
another exciting new field. Since the full RPA spectrum 
implies in general a huge number of states the mode sup- 
pression provides us a new tool for a better identification 
of distingished modes. This seems to be mandatory for 
a more reliable interpretation of fine structures in the re- 
sponse functions in comparison with measured data. 
The author acknowledges thankfully worthwile discus- 
sions with Dr. R. G. Nazmitdinov. 
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